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ABSTRACT 

Context. Bright-rimmed clouds (BRCs) are potential examples of triggered star formation regions, in which photoioni- 
sation driven shocks caused by the expansion of HIT regions induce protostellar collapse within the clouds. 
Aims. The main purpose of the paper is to establish the level of star formation occuring within a known set of BRCs. 
A secondary aim is to determine the extent, if any, to which this star formation has been promulgated by the process 
of photoionisation triggering. 

Methods. A primary set of observations is presented obtained with submillimeter SCUBA observations and archival 
data from near-IR and mid- to far-IR have been explored for relevant observations and incorporated where appropriate. 
Results. SCUBA observations show a total of 47 dense cores within the heads of 44 observed BRCs drawn from a 
catalogue of IRAS sources embedded within HII regions, supportive of the scenario proposed by RDI models. The 
physical properties of these cores indicate star formation across the majority of our sample. This star formation ap- 
pears to be predominately in the regime of intermediate to high mass and may indicate the formation of clusters. IR 
observations indicate the association of early star forming sources with our sample. A fundamental difference appears 
to exist between different morphological types of BRC, which may indicate a different evolutionary pathway toward 
star formation in the different types of BRC. 

Conclusions. Bright-rimmed clouds are found to harbour star formation in its early stages. Different evolutionary sce- 
narios are found to exist for different morphological types of BRC. The morphology of a BRC is described as type 
'A', luoderately curved rims, type 'B', tightly curved rims, and 'C, cometary rims. 'B' and 'C morphological types 
show a clear link between their associated star formation and the strength of the ionisation field within which they are 
embedded. 

An analysis of the mass function of potentially induced star-forming regions indicate that radiatively-driven implosion 
of molecular clouds may contribute significantly toward the intermediate to high-mass stellar mass function. 
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1. Introduction protostars associated with HII regions are systematically 

more luminous than those that are not ([Dobashi et aLl 

Bright-rimmed clouds (BRCs) are small molecular clouds 2OOI,). This highlights a trend toward higher mass star for- 

found at the edges of large HII regions. These clouds nation and greater star formation efficiency. Confirming 

are potential examples of triggered star-forming regions, bright-rimmed clouds as star-forming and quantifying the 

whereby shocks driven into the BRC by the photoiomsa- ^^^^^^ ^^^ location of their star formation thus provides 

tion of the cloud exterior results m the collapse (and per- important insights into the clustered mode of star forma- 

haps even the formation) of sub-critical molecular cores ^^^^ ^nd the overall star formation efficiencies of molecular 

within the cloud (e .g. Dehar veng et al.l |2005; Elmcgreen! clouds 
1X9981 : ISugitani et aLiriQQI ) . The process of photoionisation- 

induced coll apse is usually known as radiative-driven implo - ISugitani etai\ (|l99ll ). hereafter referred to as SF091, 

sion or RDI (|Bertoldilll989l:lLefloch fc Lazareflll995lll994l^. searched the Sharpless HH region catalogue (ISharolessI 

The RDI of molecular cores at the periphery of HII regions 1959) for bright-rimmed clouds associated with IRAS 

may thus be responsible for subsequent generation of star point sources, their selection was based upon the FIR 

formation, amounting to a possible cum ulative total of sev - colours thought to id entify YSOs/protostars. Later, 

eral hundred new stars per HII region (jOgura et a l.l l2002[ ) ISugitani fc Qgural (|1994D - S094 - extended their search 

and perhaps 15% or r nore of the low- to-intermediate mass to include bright-rimmed clouds from the ESO(R) Southern 

stellar mass function (jSugitani et al.l[T991.) . IRAS-selected Hemisphere Atlas. A total of 89 optically identified bright- 

rimmed clouds have been found to be associated with 

Send offprint requests to: L.K.Morgan IRAS point sources. For brevity (and consistency with 



L.K.Morgan et al.: A SCUBA survey of bright-rimmed clouds 



SIMBAD), we will refer to the combined SF091 and S094 
catalogues as the SFO catalogue. Whilst a few individ- 
ual clouds from the SFO catalogue have been studied 
in detail fe.g. iLefloch et all llQQTt iMeeeath fc Wilsonl 



1997t ICodella et al.l 120011 [Thompson et all l2004b[ 
Thompson fc Whitd 120041 : lUrguhart et all 12001 12006 . 
2007| ) and shown to harbour protostellar cores, the ques 



tion of whether star formation is a common occurence 
within BRCs (irrespective of the formation mechanism) is 
still unresolved. 

We are in the process of carrying out a star formation 
census of SFO bright-rimmed clouds in order to determine 
the current level of star formation within these sources 
and to relate it to their physical properties and morpholo- 
gies with respect to their ionising star(s). In two previous 
publications we reported observations o f the ionised gas 
surrounding the bri ght-rimmed clouds (JThompson et ahl 
l2004at [Morgan et al.ll2004 ) . Here we report a Submillimetre 
Common User Array (SCUBA) imaging survey of the sub- 
millimetre continuum emission from the clouds, in order 
to reveal the presence of compact, potentially star-forming, 
dust cores within the clouds and to contrast the properties 
of these cores with those found in other star-forming re- 
gions. We use NIR archival data from the 2 millimetre all- 
sky survey (2MASS) database to construct colour-colour 
diagrams to investigate the current level of star formation 
within these clouds and to identify YSOs. We examine the 
kinematics and temperature of the dust cores revealed by 
our SCUBA observations to show that a number of the 
SFO bright-rimmed clouds are associated with active star 
formation. 



2. Observations 

2.1. SCUBA continuum maps 

Simultan eous 450 and 850 /xm images were obtained using 
SCUBA (jHolland et al.lll999f) on the James Clerk MaxweU 
telescope ( JCMlQ) . SCUBA is a dual-camera system com- 
prised of two bolometer arrays which simultaneously sam- 
ple a similar field of view (~ 2' square). Each bolometer 
array is arranged in a hexagonal pattern and both arrays 
are observed using a dichroic beamsplitter. The 91 pixel 
short-wave array is optimised for operation at 450 /zm and 
the 37 pixel long-wave array of 37 pixels is optimised for 
operation at 850 /zm . The arrays do not fully spatially sam- 
ple the field-of-view and in order to provide fully-sampled 
images the secondary mirror of the JCMT is moved in a 
64-point pattern ("jiggling") whilst also chopping at a fre- 
quency of 1 Hz to remove sky emission. This procedure is 
commonly known as a "jiggle-map" . 

Our initial source list was formed from those BRCs in 
the SFO catalogue lying within the declination limits visible 
from the JCMT, which gives a total number of 50 clouds (45 
from the mainly northern hemisphere SF091 catalogue and 
5 from the southern hemisphere S094 catalogue). We also 
included the two bright-rimmed clouds SFO IINE and SFO 
HE in our sample. These two clouds li e in close proximit y 
to the bright-rimmed cloud SFO 11 (jOguraet al.ll2002D . 
SCUBA observations of these three clouds have already 



^ The JCMT is operated by the Joint Astronomy Centre 
on behalf of PPARC for the United Kingdom, the Netherlands 
Organisation of Scientific Research and the National Research 
Council of Canada. 



been reported in [Thompson et all (|2004bl ). Thus our initial 
source list comprised 52 bright-rimmed clouds. However, 
observations of only 47 clouds were completed, due to the 
fact that our observations were performed in the flexibly- 
scheduled mode used at the JCMT. As observations of 
three of these c l ouds h ave been previously reported by 
iThompson et al.l (|2004bf ) we present here observations of 
the remaining 44 clouds. Flexibly scheduled observations 
are not scheduled over a pre-defined period but are car- 
ried out according to the appropriate weather (atmospheric 
opacity) band, the visibility of the sources from the tele- 
scope and the scientific priority of the observations. The 
observations of the BRCs in our sample took place over 
several days from November 2001 to April 2002. 

As the angular diameters of most of the clouds in the 
SFO catalogue are comparable to the ficld-of-view (FOV) 
of SCUBA we obtained single jiggle-maps for each cloud. A 
small number of clouds in the sample either have larger 
angular diameters or lie in associations that are larger 
than the SCUBA FOV. In these cases multiple overlap- 
ping jiggle-maps were obtained and mosaiced together to 
cover the cloud or association. Varying integration times 
were used for each cloud, depending upon the atmospheric 
opacity and the brightness of the submillimetre emission 
that was detected. A chop throw of 120" was initially used 
for each jiggle-map but after inspection of the first images 
this was changed to 180" to avoid chopping onto extended 
emission. The atmospheric opacity was measured using the 
CSO 225 GHz radiometer and by performing hourly skydips 
with SCUBA at the approximate azimuth of each observa- 
tion. Typical sky opacities were 2.270 at 450 fim and 0.326 
at 850 /Ltm. 

The data were reduced and calibrated in the standard 
mann er using the SCUBA User Reduction Facility SURF 
(jJenness fc Lightfootll2000f) and the Star link image analysis 
package KAPPA (jCurrie fc Berrvl [2002). Absolute flux 
calibration was performed using calibration maps of the 
primary calibrators, Uranus and Mars, or the secondary 
calibrators CRL 618, CRL 12688, 16293-2422, IRC 10216 
and OH 231.8 depending upon the observation date, time 
and sky position. Predicted fluxes for Uranus and Mars 
were estimated usin g the values given by the Starlink 
package FLUXES (jPrivett et all Il998[ ) and predicted 
fluxes for the secondary calibrators were taken from the 
JCMT calibrator webpagcQ. 

The JCMT has a well known error beam contribution 
that may be modelled by a combination of two circularly 
symmetric gaussian functions representing the main and er- 
ror beams. We removed the error-beam contribution to the 
flux by deconvolving the calibrated images with a model 
two-component beam determined from azimuthally aver- 
aged maps of the primary calibrators ( u sing the method 
outline d in iHogerheijde fc Sandelll 120001 : IThompson et al.l 
l2004b[ ). It is well known that CLEAN processes do poorly 
with extended structure and therefore as many of our 
sources are extended on the scale of the JCMT primary 
beam CLEANing was performed only until negative com- 
ponents were encountered. Images were then restored to a 
resolution of 14". This is the native resolution at 850 /xm 
and hence smooths the 450 /im data from its native resolu- 



^ http://www.jach.hawaii.edu/JACpublic/JCMT/Continuum 
_observing/SCUB A/astronomy/calibration 
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tion of 8". The smoothing has the advantages of facihtat- 
ing the comparison of the 450 /im and 850 /im images and 
improving the signal-to-noise ratio of the 450 fiia data to 
extended emission by approximately a factor of 3. 

2.2. Archival data 

2.2.1. IRAS HiRes Data 

IRAS HiRes images at 12, 25, 60 and 100 /im were ob- 
tained from the NASA/IPAC Infrared Science Archive0 to 
complement the SCUBA data, extend the FIR wavelength 
coverage of each source and enable the spectral energy dis- 
tribution (SED) to be measured. The IRAS HiRes image 
construction alg orithm utilises the m aximum correlation 
method (MCM) (lAumann et al.lll990l ) in the enhancement 
of the original all-sky survey data. The angular resolution 
and signal-to-noise ratio of HiRes data varies depending 
upon the number of deconvolution iterations used and the 
position within the image. HiRes images are also subject 
to a number of processing artefacts, most notable of which 
are the negative 'bowls' surrounding bright sources that 
are known as ringing. This effect is much more prevalent at 
higher wavelengths, we found that 60 and 100 fim images 
may need significant correction, while 12 and 25 fim images 
needed very little. Corrections to source fluxes were esti- 
mated through measuring any negative offset present close 
to the source and adjusting the source flux accordingly. 

For IRAS data sets with typical coverages, correction 
factor variance (C FV) values reach 0.01-0 .001 at strong, 
converged sources (jLevine &: Suracd I1993D . It was found 
that the default processing parameters for HiRes images 
were not sufficient to achieve this in the majority of our 
images. Visual inspection of the produced maps allowed us 
to discount the possibility that the regions were saturated, 
overly noisy or contained 'bad' scans. Coverage maps of all 
of our sources allowed us to ascertain the level of cover- 
age per source, as well as the uniformity of the coverage 
over the beam area. This was good for all sources and so 
we iterated all sources until we reached the recommended 
CFV values of 0.01-0.001, this was generally found to be af- 
ter ~150-200 iterations, similar to the number of iterations 
used in the a nalysis of the Serpens star-forming cloud core 
(|Hurt fc Bars onv 199(| 

Typical (mean) angular resolutions were 36", 34", 
49" and 73" at 12, 25, 60 and 100 /xm respectively, which 
are generally only sufficient to identify the strongest cores 
in each SCUBA field as point sources. The absolute cali- 
bration uncertainties of the HiRes data are estimated to be 
around 20%, which is similar to the original IRAS survey 
data. 

Individual fluxes for each of the SCUBA cores were mea- 
sured from the HiRes images by aperture photometry and 
are listed in Table [TJ It should be noted that the fluxes tab- 
ulated here differ from those given for each source in the 
IRAS point source catalogue. This is, no doubt, due to the 
nearly five-fold enhanced resolution of the HiRes images 
over the original IRAS survey. HiRes fluxes are typically 
2-3 times lower than those listed in the IRAS PSC with the 
largest differences seen in the 12 and 25 /Ltm bands. A com- 
parison of the luminosities as determined from HiRes fluxes 
to those derived by SF091 will be addressed in Section [531 



This is likely to have a large effect upon the conclusions 
of SF091 who concluded that the IRAS sources embedded 
within these bright-rims had large luminosities in compar- 
ison to IRAS sources embedded within dark globules and 
other dense cores. With our improved IR resolution and ex- 
tended coverage of the sources SEDs we aim to investigate 
the true luminosity function of the Sugitani catalogue. 

2.2.2. 2MASS 

Data from the nea r-infrared J, H and Kg 2MASS catalogue 
(jCutri et al.ll2003| ) were obtained to search for protostars 
and embedded young stellar objects (YSOs) associated with 
the BRCs. The catalogue data were obtained from the 
2MASS Point Source catalogue held at the NASA/IPAC 
Infrared Science ArchiveQ. The accuracy of the photomet- 
ric measurements in the 2MASS Point Source Catalogue 
is between 1 and 2%, though can be much larger for weak 
sources. The limiting magnitudes at J, H and Ks in the 
2MASS catalogue are 15.8, 15.1 and 14.3 respectively. 

3. Results 8i Analysis 

Of the forty-four BRCs that were observed, 39 yielded a 
detection of at least one submillimetre core and five showed 
no evidence of submillimetre emission. We define a core as a 
closed contour of submillimetre emission at or above the 3a 
level and equal to or greater than the size of the FWHM of 
the JCMT beam. Individual cores were identified through 
a visual inspection of the submillimetre images and while 
most of the BRCs were found to be associated with a single 
submillimetre core, a small number of clouds were found to 
possess multiple cores (SFO 25, SFO 39 and SFO 87). 

Clouds that did not exhibit any detectable submillime- 
tre emission greater than 3 times the r.m.s. flux level were 
automatically classified as non-detections and all sources 
were assessed visually to determine if the map showed emis- 
sion consistent with an embedded core and if that emission 
was spatially correlated between the two wavelengths. If the 
source failed these two criteria then it was also classified as 
a non-detection. At 850 ^m SFO 06, 17, 28, 77 and 78 failed 
these criteria and 19 clouds failed these criteria at 450 /im. 
Upper limits to the fluxes are given in Tabled 

We present the SCUBA images as contour plots overlaid 
on DSS R images in the online appendix. The properties of 
the individual cores are described in Section [3.1l 

3.1. Core Positions, Fluxes and Sizes 

The position of each core, their measured fluxes and effec- 
tive diameters are given in Table [2] The positions given for 
each core are those determined from the peak of the 850 
/im emission. Multiple cores have been designated with an 

SMM number to indicate each compo nent. 

The Starlink package GAIA (Dr aper et all |2004[) was 
used to measure the peak and integrated fluxes of each 
core, where integrated fluxes represent a summation over 
a 30" diameter aperture centred at the position of peak 
flux, a 30" aperture was chosen as it represents a width 
of 5(7 of a 14" beam. Background levels were estimated 
from emission- free regions of each map and subtracted from 
the measured flux values. We estimate that the systematic 



http://irsa.ipac.caltech.edu 
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Table 1. Fluxes of each source as measured in IRAS HiRes images. 



Source 




IRAS HiRes Flux(Jy) 




12 ^m 


25 fiva 


60 fiva 


100 /xm 


SFOOl 


1.3 


2.9 


46.8 


173.2 


SFO 02 


0.9 


1.5 


13.6 


51.9 


SFO03 


0.2 


1.5 


12.2 


52.1 


SFO 04 


0.5 


0.9 


11.9 


25.2 


SFO 05 


1.8 


30.2 


95.0 


155.4 


SFO 06 


0.2 


0.3 


3.4 


15.5 


SFO 07 


0.5 


1.1 


18.0 


74.9 


SFO 09 


0.4 


0.6 


15.5 


47.6 


SFO 10 


1.0 


1.9 


25.6 


72.8 


SFO 12 


0.4 


1.1 


16.7 


36.6 


SFO 13 


1.2 


2.3 


35.9 


84.5 


SFO 14 


5.2 


49.6 


420.2 


534.7 


SFO 15 


0.2 


1.0 


7.3 


23.6 


SFO 16 


0.1 


1.3 


14.6 


23.0 


SFO 17 


0.1 


0.7 


3.3 


16.0 


SFO 18 


0.1 


1.0 


14.0 


22.2 


SFO 23 


0.1 


0.4 


4.1 


9.4 


SFO 24 


0.3 


0.6 


8.7 


26.0 


SFO 25 


0.3 


1.4 


17.9 


44.1 


SFO 26 


0.2 


0.3 


4.7 


17.7 


SFO 27 


0.5 


1.8 


15.8 


97.2 


SFO 28 


1.4 


1.6 


37.4 


126.6 


SFO 29 


0.3 


0.3 


5.7 


29.5 


SFO 30 


9.1 


69.0 


665.0 


901.2 


SFO 31 


1.0 


3.8 


26.2 


127.5 


SFO 32 


0.4 


0.2 


6.8 


25.7 


SFO 33 


0.2 


0.8 


7.0 


43.4 


SFO 34 


0.6 


1.0 


6.8 


24.7 


SFO 35 


0.3 


0.6 


4.2 


22.7 


SFO 36 


0.3 


1.7 


14.9 


54.8 


SFO 37 


1.8 


11.2 


35.6 


45.2 


SFO 38 


0.9 


3.6 


62.2 


172.2 


SFO 39 


0.7 


3.8 


17.4 


33.7 


SFO 40 


0.5 


0.4 


5.6 


25.8 


SFO 41 


0.5 


0.7 


5.7 


21.6 


SFO 42 


0.3 


0.6 


5.6 


20.1 


SFO 43 


2.8 


12.1 


95.4 


144.0 


SFO 44 


0.6 


8.7 


124.4 


300.2 


SFO 45 


0.1 


0.3 


1.2 


7.7 


SFO 77 


0.5 


5.1 


55.7 


67.1 


SFO 78 


0.4 


1.1 


9.7 


37.2 


SFO 87 


3.3 


4.2 


138.5 


430.6 


SFO 88 


1.6 


3.8 


120.7 


396.4 


SFO 89 


3.8 


5.5 


111.9 


330.2 



errors in measuring the fluxes of the cloud cores are no 
more than 30% in the case of the 450 /im measurements 
and 10% for the 850 /im measurements (including errors in 
the absolute flux calibration). 

Because of the higher signal-to-noise ratio of the 850 
fiin maps in comparison to those at 450 fim the effective 
angular diameter of each source was estimated from the 
longer wavelength maps. A Gaussian function was fitted 
to the azimuthally averaged 850 fim flux. The beam size 
was taken into account through the assumption of a simple 
Gaussian source convolved with a Gaussian beam of 14" 
(0obs=0bcam+0?ourcc); Three sourccs (SFO 33, 41 and 45) 
were found to be marginally resolved or unresolved at 850 
/im . 

The distances to each cloud given in SF091 and S094 
allowed a calculation of the physical effective diameter of 
each core (see Table [2]). These spatial diameters range from 
0.01-0.71 pc with a mean of 0.17 pc. 



The size of the cores may be placed in the context of star 
formation by considering turbulent motions within molecu- 
lar clouds. Observed turbulent motions become subsonic on 
the smallest scales creating smooth, regular structures. The 
size of the smallest star-forming cores may then be deter- 
mined by the scale at which turbulence transits from super- 
sonic to subsonic. This is of the order 0.05 - 0.1 pc (jLarsonl 
2003) (the resolution of the JCMT, 14", corresponds to 
0.09 pc at the mean distance of the sample, 1.29 kpc). The 
cores in this sample are evenly distributed in terms of size 
around this turbulent transitional scale. The correlation of 
the core sizes with the transitional size scale is supportive 
of the identification of the cores as being star-forming or, 
having the potential to be such. 
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Table 2. Core positions, fluxes and effective diameters. Upper peak flux limits for non-detections are given as the three- 
sigma value. The integrated detection limit is taken as the three-sigma flux value integrated over our 30" aperture. Where 
there are no detections at either 450 or 850 /Ltm the superscript " indicates that the recorded position is that of the map 
centre. 













Integrated Flux 








12000 


(^2000 


Peak Flux(Jy/beam) 


(Jy) 


Dett 


Rim Type 


Source 


(Peak) 


(Peak) 


450 fim 


850 p,m 


450 Aim 850 /^m 


(pc) 


fSeetionl4.4t 


SFO 01 


23:59:32.3 


+67:24:03 


2.8+0.2 


0.69+0.07 


8.4+2.3 0.99+0.19 


0.03 


B 


SFO 02 


00:03:58.6 


+68:35:12 


2.8+0.2 


0.60+0.03 


15.6+2.0 1.03+0.07 


0.10 


A 


SFO 03 


00:05:22.7 


+67:17:56 


<1.5 


0.29+0.03 


<13.7 0.47+0.09 


0.08 





SFO 04 


00:59:01.0 


+60:53:29 


2.7+0.4 


0.16+0.02 


10.3+3.0 0.33+0.06 


0.06 


B 


SFO 05 


02:29:02.2 


+61:33:33 


7.8+2.0 


1.44+0.04 


13.1 + 18.2 2.04+0.11 


0.16 


B 


SFO 06 


02:34:45.1 


+60:47:48 


<5.4 


<0.21 


<49.1 <0.57 


a 


A 


SFO 07 


02:34:47.8 


+61:46:29 


6.6+0.4 


1.24+0.05 


27.4+3.3 1.77+0.14 


0.19 


B 


SFO 09 


02:36:27.6 


+61:24:02 


2.4+0.4 


0.18+0.02 


6.6+3.1 0.30+0.05 


0.33 


A 


SFO 10 


02:48:12.3 


+60:24:32 


<2.6 


0.16+0.02 


<19.3 0.21+0.06 


0.35 


A 


SFO 12 


02:55:01.5 


+60:35:43 


3.1+0.7 


0.68+0.02 


1.9+4.9 1.10+0.06 


0.31 


B 


SFO 13 


03:00:56.5 


+60:40:24 


11.6+3.0 0.69+0.03 


15.6+22.3 1.33+0.09 


0.40 


B 


SFO 14 


03:01:31.3 


+60:29:20 


12.1 + 1.2 2.68+0.04 


50.8+9.0 4.54+0.11 


0.40 


A 


SFO 15 


05:23:28.3 


+33:11:48 


<1.5 


0.10+0.01 


<11.3 0.12+0.03 


0.71 


B 


SFO 16 


05:19:48.4 


-05:52:04 


±1.2 0.85+0.03 


23.0+9.2 1.54+0.07 


0.06 


A 


SFO 17 


05:31:28.1 


+ 12:05:24 


<0.7 


<0.20 


<5.4 <0.52 


a 


A 


SFO 18 


05:44:29.7 


+09:08:55 


2.0+0.2 


0.72+0.02 


6.3+1.5 1.25+0.05 


0.09 


A 


SFO 23 


06:22:57.8 


+23:10:15 


2.5+0.5 


0.23+0.02 


3.8+3.7 0.40+0.06 


0.27 


A 


SFO 24 


06:34:52.2 


+04:25:38 


<3.3 


0.19+0.02 


<21.8 0.12+0.04 


0.25 


B 


SFO 25 SMMl 


06:41:03.2 


+ 10:15:09 


5.1+0.8 


1.16+0.03 


10.0+5.2 1.66+0.08 


0.06 


B 


SFO 25 SMM2 


06:41:04.8 


+ 10:15:01 


3.7+0.8 


0.90+0.03 


10.2+5.2 1.70+0.08 


0.06 


B 


SFO 26 


07:03:46.7 


-11:45:52 


<1.4 


0.14+0.02 


<9.1 0.23+0.05 


0.13 


A 


SFO 27 


07:03:58.3 


-11:23:04 


<1.5 


0.32+0.02 


<9.8 0.59+0.06 


0.48 


A 


SFO 28 


07:04:44.4 


-10:21:41 


<1.6 


<0.14 


<15.7 <0.39 


a 


A 


SFO 29 


07:04:52.4 


-12:09:43 


<2.7 


0.45+0.06 


<25.3 0.60+0.18 


0.04 


A 


SFO 30 


18:18:46.8 


-13:44:28 


5.0+0.3 


1.89+0.03 


20.0+2.5 2.42+0.07 


0.33 


B 


SFO 31 


20:50:43.1 


+44:21:56 


2.6+0.2 


0.62+0.04 


5.6+2.3 0.72+0.12 


0.04 


A 


SFO 32 


21:32:29.5 


+57:24:33 


<0.8 


0.22+0.01 


<7.9 0.13+0.03 


0.20 


A 


SFO 33 


21:33:11.9 


+57:30:07 


<1.0 


0.14+0.02 


<8.2 0.10+0.06 


0.01 


A 


SFO 34 


21:33:32.1 


+58:03:32 


2.0+0.1 


0.48+0.01 


5.7+0.6 0.60+0.03 


0.04 


A 


SFO 35 


21:36:03.4 


+58:31:25 


<1.6 


0.20+0.02 


<10.3 0.36+0.05 


0.11 


A 


SFO 36 


21:36:07.4 


+57:26:41 


14.9+0.3 1.73+0.03 


38.7+1.7 2.12+0.07 


0.07 


A 


SFO 37 


21:40:28.9 


+56:35:53 


<1.5 


0.98+0.04 


<7.7 1.40+0.10 


0.06 


C 


SFO 38 


21:40:41.8 


+58:16:14 


6.2+0.3 


4.02+0.05 


24.0+2.0 6.68+0.14 


0.11 


B 


SFO 39 SMMl 


21:46:01.2 


+57:27:42 


3.8+0.6 


0.66+0.03 


13.5+4.0 1.12+0.09 


0.06 


B 


SFO 39 SMM2 


21:46:06.9 


+57:26:36 


3.4+0.6 


0.64+0.03 


8.7+4.0 0.94+0.09 


0.03 


B 


SFO 40 


21:46:13.4 


+57:09:42 


3.4+0.6 


0.30+0.04 


9.4+4.4 0.44+0.10 


0.10 


A 


SFO 41 


21:46:28.6 


+57:18:35 


<1.1 


0.20+0.02 


<10.4 0.38+0.05 


0.04 


B 


SFO 42 


21:46:34.7 


+57:12:31 


<1.3 


0.25+0.03 


<12.0 0.50+0.09 


0.09 


A 


SFO 43 


22:47:49.3 


+58:02:51 


2.0+0.5 


0.60+0.06 


7.6+4.4 1.19+0.16 


0.31 


B 


SFO 44 


22:28:51.1 


+64:13:40 


18.7+0.6 2.93+0.06 


63.2+5.4 3.88+0.16 


0.06 


A 


SFO 45 


07:18:26.6 


-22:05:56 


<1.1 


0.14+0.01 


<8.8 0.04+0.04 


0.05 


A 


SFO 77 


16:19:53.9 


-25:33:39 


<4.1 


<0.09 


<40.6 <0.24 


a 


A 


SFO 78 


16:20:52.9 


-25:08:07 


<3.6 


<0.10 


<35.7 <0.28 


a 


A 


SFO 87 SMMl 


18:02:49.8 


-24:22:26 


3.2+0.7 


0.64+0.04 


8.0+7.0 1.22+0.11 


0.12 


B 


SFO 87 SMM2 


18:02:27.3 


-24:22:57 


2.2+0.6 


0.44+0.02 


5.3+6.2 0.79+0.07 


0.15 


B 


SFO 88 


18:04:11.5 


-24:06:41 


1.6+0.2 


0.90+0.03 


7.8+2.3 1.96+0.09 


0.23 


A 


SFO 89 


18:09:56.6 


-24:04:21 


2.5+0.3 


0.52+0.02 


14.2+2.6 0.97+0.05 


0.21 


A 



3.2. Temperature, Mass and Density 



It has bee n shown (e.g., ISridharan et al.l l2002t 

iBeuther et al] |2002[ ) that the SED of luminous YSOs 
can be modelled by a two component modified black- 
body fit. The two components consist of a compact, hot 
component which dominates the IRAS 12 /xm and 25 
/xm fluxes and a component which represents the cooler, 
more extended, dust. The cooler component of the SED 
dominates the IRAS 60 fim and 100 /im fluxes in addi- 
tion to fluxes at submillimetre wavelengths. The cooler 
component arises from re-radiated emission, absorbed and 
then emitted by the outer envelope of the protostellar 
source. The temperature that is then associated with this 
component is referred to as the 'dust' temperature of the 
cloud. Measured fluxes classifled as non-detections have 
been discarded when fitting SEDs, thus some SEDs have 
been created using only the IRAS fluxes and the 850 
/Km flux. 

Certain sources were not subjected to SED fltting anal- 
ysis, even if they have been classified as detections. This 



was due to high levels of error in flux measurements or 
poor morphological agreement in the source structure at 
different wavelengths. These sources were SFO 24, 32, 35, 
36, 38, 44, 45 and 88. 

Fits to the IRAS & SCUBA fluxes w ere made using 
the approach outlined in (jDent et al.lll998l ) , specifically, the 
integrated flux at each frequency (F^) was determined from 
the following equation: 



F,^nB,{Td)il 



(1) 



where fl is the solid angle subtended by the aperture 
used when summing fluxes, B^(Td) is the Planck function 
evaluated at a dust temperature, Td, and frequency, v, and 
Ti, is the optical depth at frequency v. The optical depth of 

the dust was assumed to follow the form t^^ — Tj-d I -^ ) , 

where v^d is a reference frequency and /3 is the dust emis- 
sivity. As the SEDs of the cores are defined by only five or 
six points we chose to reduce the number of free parame- 
ters in our model by fixing /3. For the small, hot compo- 
nent f3 was fixed at a value of 1 in following with previous 
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works fe.g lFaundez et al.ll2004l : ISridharan et al.ll2002f) and 
at i ts canonical va lue of 2 for the cooler, more extended gas 
fe.g lThompson et al.ii2004b) . An example SED is shown in 
FigUl 

The mass of each core was determined using the method 
of lHildebrandl ([lOSS"), here we measure the mass of the cold 
components of our SEDs, the dominant component of mass 
in the cloud. For an optically thin cloud with a uniform 
temperature the total (Mdust+A^gas) mass of the cloud M 
is given by 



M 






(2) 



The parameter C^ is a mass conversion factor combining 
both the gas-to-dust ratio (typically assumed to be 100) and 
the frequency dependent dust o pacity, n^. Quoted values 
for C y range from 2 1.4 g cm~^ (IKruegel fc SiebenmorgenI 
fl99l to 286 g cm-2 (iDraine fc Ledll984D . A value of C, = 
50 g cm~^ Sit V = 850jmi has been adopted here following 
[Thompson et al.1 (|2004bl ). Number densities of the clouds 
were determined from the mass by assuming a spherical ge- 
ometry for the core. The dust temperatures resulting from 
the modified blackbody fits, along with core masses and 
H2 number densities are presented in Tabled Upper limits 
for the masses of potential cores in clouds labelled as 'non- 
detections' are presented in table [5] along with upper limits 
for density and extinction values. 

The uncertainties in the mass (and hence density) of 
the clouds are generally dominated by temperature effects 
in the non-linear Planck function in Eql^l Typical values 
of observed 850 /xm flux and SED determined dust tem- 
perature have uncertainties of ^ 10% and ± 1-2 K respec- 
tively leading to an factor of ~ v2 uncertainty in mass and 
density. In addition, it should be noted that the value of 
Cu is a factor of two smaller than that typically assumed. 
This was to maintain consistenc y with other wo rks based 
on the same set of observations (iThompson et a l. 2004b). 
This uncertainty will also propagate into our derived values 
of mass and density and so these results should be scaled 
by the best known value of C^. In addition, there will be 
some uncertainty in mass and density values introduced by 
our uncertainty in (3. This is commonly assumed to be 2 
and can be fairly well constrained in certain circumstances, 
note that only the cooler component of (3 is relevant in this 
case as we are only deriving the mass associated with this 
component of our sources, the majority of the mass of the 
system. 

Values of the visual and near-infrared extinction [Ay & 
Ak) toward each core were deri ved from their s u bmilli metre 
fluxes following the method of iMitchell et all (|200lD . The 
submillimetre flux F,j may be related to the visual extinc- 
tion by: 



Fi, = ^By{Ta)K,v'in]i 



N„ 



E{B-V)R 



(3) 



where the total opacity of gas and dust is represented 
by K^, which is functionally equivalent to the reciprocal of 
the mass conversion factor C^, mn is the molecular mass of 
interstellar material, N-h_/E{B — V) is the conversion fac- 
tor between column density of hydrogen nuclei and the 
selective absorption at the relevant wavelength. Following 



IMitchell et al.l (|2001f ) (and references therein) a value of 
Nii/E{B — V) = 5.8 X 10^^ cm~^ mag~^ has been assumed. 
A value of R = 5 has been assumed in line with predic- 
tions of th e extinction of the inne r regions of molecular 
clouds fe.g lCampeggio et al.l (120071 )1. The K-band extinc- 
tion Ak of each so urce was found by mult iplying by the 
ratio Av/Ak = 8.9 (J Rieke fc Lebofskvlll985D . Values of Ay 
and Ak are listed in Table [H along with the dust tem- 
peratures and masses, etc. derived from the submillimetre 
observations. 

3.3. 2MASS colour-selected YSOs associated with the BRCs 

The 2MASS catalogue was searched for near-infrared ob- 
jects within the optical bounds of the surveyed clouds to 
identify any YSOs that may be associated with the BRCs. 
YSOs occupy a distinctive region in a J — H and H ~ K^ 
colour diagram that is intrinsically redder than the main- 
sequence population (jLada fc Adamsl[l992D . An example of 
this type of diagram is shown in Fig. ^ with similar dia- 
grams for the rest of the clouds in the survey presented in 
the online appendix. The solid lin es represent the colours 
of main sequence and giant stars (jKoornneeJ 1 1983t ) . while 
the dotted line sho ws the Classical T-Tauri (CTTS) locus 
(jMever et al.lll997D . 

YSOs may be differentiated from older stars that have 
had their radiation extincted by large amounts of inter- 
vening dust through their infrared colours. As a source is 
subjected to higher levels of extinction (larger interven- 
ing columns of dust) it moves towards the top right cor- 
ner of the diagram, following the reddening tracks shown 
by dashed lines in Fig. [51 Thus, sources found to lie to 
the left of the reddening track associated with the latest 
main-sequence stars arc likely to be main-sequence or giant 
stars, while those sources that lie between this track and 
the rightmost track associated with the CTTS locus and 
above the locus itself may be identified as CTTS. Sources 
to the right of the rightmost track are likely to represent 
Class I protostars hidden by a large, extended dusty enve- 
lope. The numbers of CTTS and Class I protostar sources 
were determined using these constraints and are listed in 
Table H 

The regions occupied by main sequence or giant stars, 
CTTS Class II protostars and class I protostars in these di- 
agrams are not compl etely discrete; ea.ch po pulation over- 
laps with the others (|Lada fc Adamslll992D . Nevertheless 
the colour-colour diagrams indicate the general population 
of recently formed stars and protostars associated with the 
clouds and allow us to identify candidate CTTS and class I 
protostars for later followup and confirmation. When iden- 
tifying potential CTTS sources we have discounted sources 
that may simply be extincted main-sequence objects, thus, 
the number of identified CTTS sources is probably an un- 
derestimate. 

Investigation of the position of 2MASS sources in the 
plane of each DSS image fails to reveal any consistent trend 
of source position with cloud morphology, i.e. the sources 
are not consistently clustered towards either the submil- 
limetre cores or the optically bright rims. Data from the 
Spitzer telescope will yield much information upon the IR 
content of each cloud and will be presented in future work. 

We indicate the numbers of candidate CTTS and class 
I protostars associated with each BRC in Table |3l The ma- 
jority of BRCs are associated with at least one YSO and 
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Fig. 1. SED fit to ffuxes observed toward ttie object SFO f2 



in many cases multiple YSOs; this is consistent with recent 
and/or ongoing star formation activity: out of our sample 
of 47 detected cores there are 10 that do not show any 
candidate CTTS or class I protostars (SFO 1, 24, 26, 28, 
34, 78, 87a, 87b, 88 and 89). Protostellar sources in the 
early stages of evolution are likely to be highly embedded 
within their host clouds behind large degrees of extinction. 
These sources are also likely to have low fluxes at NIR wave- 
lengths, thus the number of associated 2MASS sources may 
be considered as a lower limit. This indicates that the ma- 
jority of observed BRCs are still in an early phase of their 
star formation history. The positions of each source identi- 
fied as either CTTS or protostars are overlaid on the images 
in the online appendix. The 2MASS sources do not appear 
consistently related to either the position of peak submil- 
limeter flux or the position of the optically bright rim. 



As mentioned in Section [2. 2. 2i the limiting magnitudes 
of the 2MASS are 15.8 in the J band, 15.1 at the H band 
frequency and 14.3 in the K^ band. The average if-band 
extinction of the sample set ranges from 0.2 to 4.5 with 
a mean of 1.3 (the determination of K-hawA extinction is 
discussed in Section 13. 2p . The maximum observable Kg- 
band magnitude should then be 9.8, equivalent to a 3.2 M0 
protostar with an age of 0.3 Myr a t the mean distance of 
the sample (1.29 kpc) (jZinnecker et al. 1993). The numbers 
of identified YSOs found should then accurately trace the 
star formation within the clouds, at least for protostars 
more massive than ~3.2 Mq. 

3.4. Additional IR Information 

In addition to the IRAS data presented in our SED 
analyses and the 2MASS fluxes utilised in Section 13.31 
infrared observations have been performed of large areas of 
the sky by both the Midcourse Space Experiment (MSX) 
and the Spitzer telescope. Searches were performed for 
observational data relevant to our sources in both the 
MSX and the publically available Galactic Legacy Infrared 
midplane Survey Extraordinaire (GLIMPSE) catalogues. 
Only one source from our sample was found to lie within 



the field observed by GLIMPSE, which covered only the 
central two degrees in latitude of the Galaxy. This source 
was SFO 30 and will be the subject of a more detailed 
study by the authors in the future. Many of our sources 
were found to have been observed at 24 ^m in proprietary 
Spitzer programmes and these data were examined for 
evidence of protostellar content within our sample. A total 
of 37 multiband imaging photometer for Spitzer (MIPS) 
maps were obtained from the Spitzer archive, covering 
sources from the SFO catalogue. Of the sources covered 
by these maps four (SFO 8, 19, 21 and 22) were not 
observed by the reported SCUBA observations and an 
additional two (SFO 6 and 17) were non-detections in the 
SCUBA survey. Of the remaining 31 sources, 25 show 24 
/im emission clearly associated with the SCUBA cores 
identified in Sect ion [XT] The other six clouds appear to 
have non-correlated 24 /xm and 850 /im emission. 

The correlation of 24 /im emission with the positions of 
the reported 850 /im cores highly supports the identification 
of these cores as star-forming. This includes cores in which 
no 2MASS YSO sources were identified (Section l^X^ . 
these sources are SFO 01, 24, 34 and 88. The non-detection 
of these cores by the 2MASS is most likely due to these 
sources being at early stages of their evolution as described 
in Section l2.2.2l A more detailed analysis of the publicly 
available Spitzer archival data is underway by the authors 
and will be reported in a future work. 

Several of our sources were found to contain objects 
identified in the MSX survey. However, a very limited num- 
ber of these sources were found to be closely associated with 
the submillimeter cores identified in this work (the criterion 
being a position within one IRAS beam- width). An even 
smaller number of these sources were found to have reliable 
fiuxes within the MSX catalogue. Sources with flux quali- 
ties of lower than a 5-sigma detection were rejected. After 
applying these criteria, sources were found associated with 
SFO 5, SFO 31, SFO 39 and SFO 43. 

Fluxes from the C and E bands of the MSX survey 
were plotted upon the SEDs of the relevant sources (SEDs 
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Fig. 2. J-H versus H-Kg diagrams of the 2MASS sourc es associated with the cloud SFO 12. The sohd hues represent the 
unredde ned loci of main seq uence and giant stars from lKoornneeJ (jl983() The dotted line represents the classical T-Tauri 
locus of lMever et al.l (jlQQTT ). Reddening tracks are shown by dashed lines and the dash-dot line represents the reddening 
track of the earliest main-sequence stars. 

Table 3. Numbers of class I protostellar candidates (CI) or candidate T-Tauri stars (CTTS) identified from the JHKg 
diagrams (online appendix) for all sources 



Source 


Number of 


Number of 


Source 


Number of 


Number of 




associated CTTS 


associated CI 




associated CTTS 


associated CI 


01 








30 


2 


1 


02 


3 


10 


31 


1 





03 





1 


32 





2 


04 





2 


33 





2 


05 


2 


1 


34 








06" 





6 


35 





3 


07 





4 


36 





3 


09 





1 


37 





4 


10 





2 


38 


1 


3 


12 


4 


3 


39 





8 


13 


1 


3 


40 





2 


14 


1 


9 


41 





2 


15 





1 


42 





1 


16 





3 


43 


1 


2 


IT 


1 





44 


1 


2 


18 





2 


45 


1 


3 


23 





2 


77" 





1 


24 








78" 








25 


2 


3 


87a 








26 








87b 








27 


4 


2 


88 








28'' 








89 








29 


1 


1 









^Non-detections at SCUBA wavelengths 



for SFO 31 and SFO 39b are shown in FigE]). These fluxes 
(found at wavelengths of 12.13 and 21.34 /im) serve as a use- 
ful comparison to our IRAS HiRes fluxes, found in Section. 
[2XT1 The SEDs show that MSX fluxes largely agree with 
the IRAS HiRes fluxes found at similar wavelengths for 
three of the five sources found to have relevant detections 
in the MSX survey, given the errors inherent in compar- 
ing the two surveys. This similarity is exemplified by the 
source SFO 31, shown in FiglHl Attention may be drawn to 
the two sources within which fluxes were not in such good 
agreement. SFO 39a and 39b show rather large disagree- 
ment between the MSX-derived fluxes and the line of the 
SED as drawn by our best-fit method. It should be noted 
that these objects are fairly poorly fit in our SED analysis 
and that the IRAS HiRes 60 /.tm point is also discrepant 
with the line. A possible explanation for this is that the 



source is known to comprise at least two distinct sources 
and so there may be 'blending' of the source emission in at 
least the IRAS beam if not the SCUBA beam. As further 
Spitzer data is released, the quality of NIR measurements of 
these sources will allow a great improvement in the fitting 
of these sources' SEDs at shorter wavelengths. 



4. Analysis 

4.1. Cloud Morphology 

The SCUBA images of the cloud sample show submil- 
limetre cores generally located immediately behind the 
optically bright rim (online appendix). The submillimetre 
core positions are generally consistent with the IRAS 
source positions presented in SF091 and S094 (i.e. the 
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Fig. 3. SEDs of the sources SFO 31 (left) and SFO 39b (right) including fluxes from the MSX survey. MSX fluxes are 
shown as error bars only, IRAS and SCUBA fluxes are plotted as diamonds. 



submillimetre cores lie within the IRAS beam 'footprint'), 
though some discrepancy may be due to the IRAS tracing 
the more diffuse dust surrounding the cores. 

The submiUimetre emission from each core is approx- 
imately radiaUy symmetric at the 50% level. Below the 
50% flux level the weaker, more diffuse, submillimetre emis- 
sion from each source generally follows the optical bound- 
aries of the bright rimmed clouds and thus may show a 
cometary morphology (e.g. SFO 13, Fig. 2]). A general sum- 
mary of morphology is then: radially symmetric contours 
surrounded by diffuse emission following the optical cloud 
morphology and often elongated along the direction of the 
ionising source, exemplified by SFO 12 in Fig. [5] 



4.2. Luminosities of Embedded Sources 

Use of the far infrared (FIR) and submillimetre fluxes of 
the sources allows the estimation of the spectral class of any 
YSOs that may be present in the cloud cores. The luminosi- 
ties of each source were estimated by integrating under the 
modified blackbody curve fitted to each cloud's observed 
fluxes, as presented in Table [H Distances to each source 
were taken from SF091 and S094. The derived properties 
of the internal sources are presented in Table |4l upper flux 
limits from those sources that were not detected by the 
survey have been utilised to determine upper limits for the 
masses of potential cores in those clouds, along with upper 
limits for density and extinction values. These values are 
presented in Table (3] 



The assumption has been made in this analysis that all 
FIR and submillimetre luminosity is due to a single source, 
enabling an order of magnitude estim ate of the stellar con- 
tent w ithin the protostellar envelope. IWood &: Churchwelll 
(|1989( ) showed that, for a realistic mass function, the spec- 
tral type of the most massive member in a cluster is only 
1.5-2 spectral classes lower than that derived for the sin- 
gle embedded star case. The median luminosity of internal 
sources in the sample is ~63 L©, corresponding to a spec- 
tral type of B9. 



4.3. Classification of the Internal Sources 



lAndre et al.l (|199 3^ defined Class protostars as those 
obeying the approximate relationship Lsubmm/Lboi ^ 5 
X 10"'^ (here we adopted the definition of Lsubmm as the 
luminosity radiated longward of 350 yum). The clouds that 
have been modelled using the modified blackbody analysis 
all approximately obey this relationship (see Table 2]) 
with 29 of 34 sources having Lsubmm/Lboi > 5 x 10^^ 
and all sources having Lsubmm/Lboi > 2.3 x 10"'^. No 
correlation was found between the masses, luminosities 
or temperatures of the sources and the value of their 

Lsubmm/Lboi ratlO. 

The dust temperatures that have been found for the 
cores in the sample presented here are higher tha n that 
normally found for starless cores (Td ^ 10 K; () Evans! 
1999)). This is expected as internal hot components are 
observed, these are heating the protostellar cores from 
within, and this establishes a scenario of ongoing star 
formation within our sample. 



The bolometric luminosities derived for the clouds in 
this survey range from 1 L©, typi cal of class and class 1 
proto stellar objects ([Andre et al.lll993t IChandler &: Richer] 
I2OOOI) to 6945 L(7i, more indicative of massive star forma- 
tion ([Mueller et al.ll20d^ . The median luminosity found for 
the sample is ~63 L©. This suggests that the sample con- 
sists largely of intermediate-mass star forming regions or 
stellar clusters, c o nsiste nt with the findings of SF091 and 
lYamaguchi et al.1 (119991 ). This is despite the differences in 
fluxes found between the IRAS PSC fluxes used by SF091 
and IRAS HiRes fluxes used here. 



4.4. Cloud Classification 

SF091 separated their catalogue into three morphologi- 
cal types based upon the curvature exhibited by the BRC. 
Their classifications are as follows, (1) Type A, moderately 
curved rim with a length to width ratio, l/w, of less than 
0.5; (2) Type B, tightly curved rim with l/w greater than 
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Table 4. Core properties determined from modified blackbody fits. 





















Distance 




Source 


Tnot 


Tcold 


M 


Logio (n) 


Av 


Ak 


Luminosity 


Spectral 


To Cloud'' 


-'-'subniiTi/-'-'bol 




(K) 


(K) 


(Mq) 


cm~^ 






(L©) 


Type" 


(kpc) 


(x 10-3) 


01 


169 


26 


2.18 


4.6 


8.5 


1.0 


122 


B8 


0.85 


6.9 


02 


182 


22 


3.37 


4.8 


13.2 


1.5 


57 


B9 


0.85 


15.5 


03 


126 


24 


1.15 


4.3 


4.5 


0.5 


39 


AO 


0.85 


9.2 


04 


182 


26 


0.03 


4.6 


2.1 


0.2 


1 


G2 


0.19 


6.7 


05 


107 


24 


25.85 


4.6 


20.2 


2.3 


1154 


B4 


1.90 


6.5 


07 


170 


21 


31.26 


4.7 


24.4 


2.7 


400 


B5 


1.90 


19.0 


09 


185 


27 


3.03 


3.7 


2.4 


0.3 


181 


B7 


1.90 


6.5 


10 


181 


29 


2.22 


3.5 


1.7 


0.2 


268 


B6 


1.90 


4.0 


12 


157 


21 


16.87 


4.4 


13.2 


1.5 


201 


B7 


1.90 


18.5 


13 


177 


23 


19.60 


4.5 


15.3 


1.7 


444 


B5 


1.90 


11.5 


14 


121 


27 


51.79 


4.9 


40.5 


4.5 


3205 


B2 


1.90 


5.8 


15 


130 


25 


7.04 


3.3 


1.7 


0.2 


317 


B6 


3.40 


7.3 


16 


60 


17 


1.25 


5.3 


22.1 


2.5 


8 


A8 


0.40 


24.4 


18 


108 


18 


1.15 


5.3 


20.3 


2.3 


6 


FO 


0.40 


30.6 


23 


125 


19 


6.26 


4.2 


6.9 


0.8 


40 


AO 


1.60 


27.1 


25a 


139 


19 


5.59 


5.1 


25.9 


2.9 


43 


AO 


0.78 


24.3 


25b 


139 


19 


5.69 


5.1 


26.4 


3.0 


43 


AO 


0.78 


24.6 


26 


212 


24 


1.05 


3.9 


2.2 


0.3 


28 


Al 


1.15 


10.2 


27 


153 


24 


3.36 


4.4 


7.2 


0.8 


108 


B8 


1.15 


9.9 


29 


219 


22 


3.57 


4.4 


7.6 


0.9 


48 


AO 


1.15 


14.5 


30 


127 


34 


29.33 


4.5 


17.1 


1.9 


6945 


B1.5 


2.20 


2.3 


31 


147 


27 


1.19 


4.1 


3.4 


0.4 


96 


B8 


1.00 


5.3 


33 


149 


24 


0.45 


4.1 


2.2 


0.3 


21 


A2 


0.75 


10.0 


34 


188 


21 


1.55 


4.6 


7.8 


0.9 


23 


A2 


0.75 


15.7 


37 


127 


20 


3.82 


5.0 


19.2 


2.2 


63 


B9 


0.75 


12.5 


39a 


133 


20 


2.77 


4.9 


13.9 


1.6 


36 


AO 


0.75 


16.3 


39b 


134 


21 


2.09 


4.7 


10.5 


1.2 


35 


AO 


0.75 


13.1 


40 


224 


24 


0.45 


4.1 


2.2 


0.3 


16 


A4 


0.75 


9.4 


41 


189 


24 


0.43 


4.0 


2.1 


0.2 


15 


A4 


0.75 


8.5 


42 


187 


22 


0.89 


4.4 


4.5 


0.5 


16 


A4 


0.75 


12.3 


43 


143 


26 


28.34 


4.4 


13.9 


1.6 


1385 


B3 


2.40 


6.1 


87a 


204 


30 


5.93 


4.4 


8.8 


1.0 


698 


B4 


1.38 


4.0 


87b 


204 


31 


4.45 


4.3 


6.6 


0.7 


657 


B4 


1.38 


3.4 


89 


194 


30 


5.18 


4.3 


7.7 


0.9 


597 


B5 


1.38 


4.1 








'^ Spectral types taken fromlDe Jager & Nieuwenh 


uijzcn I1987D 










b 


Distances 


taken from|Sugitani 


et al.l (119911) andlSueitani & Osura| l|l994|l 





Table 5. Limiting values of potential core properties determined from upper flux limits in the survey's non-detections. 



Source 



M 
(Ms) 



Logio(n) 



A^ 



Ak 



0.5 and; (3) Type C, cometary rini^- 
The three rim types are closely related to the develop- 
ment of clouds under the influence of ionising sources. The 
three general rim types, A, B and C, closely match the 
iLefloch fc Lazarefi |l994) RDI model snapshots at 0.036, 
0.126 and 1.3 Myr respectively. This strongly suggests 
that the variation in the large scale morphology repre- 
sents the duration of each clouds exposure to ionising ra- 
diation, rather than a fundamental difference in structure 
or composition. T his is further supported by the fact that 
iLefloch fc Lazarefil ([1994) found no significant effects upon 
the morphological evolution of the clouds due to varying 
initial conditions. 



06 


7.80 


4.1 




6.1 


0.7 


17 


0.32 


4.7 




5.6 


0.6 


28 


1.95 


4.1 




4.2 


0.5 


77 


2.48 


3.8 




2.6 


0.3 


78 


2.89 


3.8 




3.0 


0.3 








A 


simplistic 



interpretation of the results of 
ILefloch fc Lazareff (1994) is that, if their evolution- 
ary scenario is correct, then the fractional number of type 
'A', 'B' and 'C type rims observed should match the frac- 
tional lifetime of each phase. This was found not to b e the 
case, however. The simulations of iMiao et al.l (120061) have 
reprod uced the morphological results of ILefloch fc Lazarefil 
but incorporate self-grav ity and thermal evolu- 



tion into their models. Whilst iMiao et all (|2006f) find 



•" SFO 36 was identified by SF091 as 'A or C. Based upon 
the position and morphology of the submillimetre emission pre- 
sented in these images the identification of SFO 36 as Type 'A' 
morphology is evident (see Fig. |6} . 



a sim ilar morphological evolution to iLefioch fc Lazaren 
()1994| ) they find that the inclusion of self-gravity de- 
creases the timescale for the cometary morphology of 
their simulated clouds to dev elop. The simulations of 
iKessel-Devnet fc BurkertI (|2003[ ) include the effects of self- 
gravity as well as ionisation. These simulations similarly 
show decreased timescales for morphological evol u tion i n 
comparison to the results of ILefloch fc Lazarefl 



They also show fragmentation occuring in BRCs once a 
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type 'C morphology has been achieved (on a timescale 
of ^^0.4 Myr) and subsequently show small 'trunk-like' 
str uctures eme rging from the larger BRC. 

iMiao et all (|2006D gives representative ages for the 'A', 
'B' and 'C type rims of 0.24, 0.30 and 0.3 9 Myr respec- 
tively. Using the results of iMiao et al.l (|2006D a close match 
was found between the number of 'A' type rims classified in 
the SFO catalogue and the lifetime of each of these phases 
(63% predicted compared to 66% observed). If the varia- 
tion in rim morphology seen across the SFO catalogue may 
then be attributed to an evolutionary sequence then it is 
reasonable to search for differences in the physical proper- 
ties of the different rim types. It is expected that the core 
properties such as temperature and density would increase 
with time, while the mass of the cold component would 
decrease as the protostellar envelope is accreted onto the 
central object and the cloud is further evaporated by the 
ionising UV radiation. One would naively expect to find 
star formation that is more developed in the clouds that 
are further evolved. 

The morphological types of each BRC as identified 
by SF091 are tabulated in table [H It can be seen that 
the number of 'A' type clouds exceeds the number of 'B' 
types by approximately a factor of 2 (26 as compared to 
14) with only 2 'C type rims. While the timescales re- 
ported for 'A' type ri ms versus the tot al of 'B' and 'C 
type rim timescales in iMiao et al.l (J2006t) match those ob- 
served, the lack of observed 'C type rims is unexplained. 
Given the large amount o f time attributed to this phase by 
iLefloch fc Lazaren (|1994| ) many of these type rims should 
be observed. However, without more detailed information 
on the modelling of the eventual fate of these clouds, it 
is not possible to draw any conclusions about the lack of 
observed 'C type rims. 

Grouping the cores in the sample according to the clas- 
sifications of SF091 (TableHj) does not show any significant 
differences in the derived core properties (e.g. mass, den- 
sity and extinction) between the 'A' rims and the 'B' rims. 
This is not what might be expected if the rim types are in- 
deed an evolutionary sequence. This suggests two possible 
scenarios. 

1 . The timescale for the development of rim type is short 
compared to star formation. Thus, assuming that the 
stars are formed as a result of the interaction of the 
ionisation front with the BRC, these stars begin to 
form quickly as a result of this influence. The clouds 
then rapidly evolve through the suggested development 
of rim type irrespective of the state of star formation 
within them. 

2. The formation of stars within the BRCs is unrelated to 
the interaction of the ionisation front with the BRC. 
In this case the development of stars within the clouds 
would be unconnected to the wider cloud evolution. 



4.5. Radio Detections and Ionising Stars 



Morgan et all (|2004[ ) used archival NVSS (|Condon et all 
1998[) data to look for radio emission from the northern 



SFO catalogue. Of the 44 clouds observed here, 39 were 
observed and 31 detected by the NVSS (All of the northern 
sources tha t have been observe d with SCUBA here were 
included in iMorgan et al.l (|2004l )). All but one of these de- 
tections was classified as an ionised boundary layer (IBL) 



at the edge of the respective bright rim, establishing the in- 
teraction of the HII region with the surrounding molecular 
material. The number of sources with a clearly defined IBL 
demonstrates that ^70% of these BRCs are being influ- 
enced by the expansion of their respective ionisation fields. 
Here we find eight clouds in which we have detected submil- 
limetre emission, associated with a prospective protostellar 
core, which show no detectable radio emission. These clouds 
are SFO 3, 9, 23, 24, 26, 33, 34 and 39. 

If an interaction between the ionising source and the 
molecular material exists in the eight clouds not detected at 
~ 21 cm, then it is presumably at a much lower level than in 
the rest of the sample. A comparison between these two sets 
of sources is therefore warranted to search for differences 
which may highlight the presence and/or effects of RDI. 
The mean temperatures, masses, densities, extinctions, lu- 
minosities, internal stellar contents and the degree of IR 
source association of the two sets were compared (Table [7|). 
Clouds without associated radio emission show a tendency 
towards lower mass, density, visual extinction and luminos- 
ity, though this tendency must be dealt with cautiously due 
to the size of the standard deviations within the respective 
groups. It can be seen from Table [7] that all values agree 
within these standard deviations. A tentative statement is 
then, the presence of a strong ionisation field in the vicin- 
ity of protostars taken from the observed sample may lead 
to more massive, luminous protostellar sources. Support for 
this statement may be found by comparing the distribution 
of source mass for those sources detected and not-detected 
at radio wavelengths. Figure [7] shows a histogram of both 
sets of sources, it can be seen that the distribution of those 
sources detected at radio wavelengths contains a noticeably 
higher number of higher mass sources. 

If there exists a correlation between the radio flux em- 
anating from a source rim and the mass/luminosity of the 
internal protostar embedded within that rim then it is war- 
ranted to compare the luminosities and masses of each 
source to the predicted ionising flux of each rim's ionising 
source, i.e. the ionising flux we expect to be arriving at each 
cloud's rim given the spectral type of the star powering the 
related HII region and the projected distance between that 
star and the observed bright rim. 

A correlation was found between the predicted ionis- 
ing flux and the luminosity of each protostellar object. 
Futhermore, it was found that this correlation appeared 
to arise solely from the type 'B' and type 'C rims within 
our sample. Type 'A' rims showed no such correlation, 
a plot of the relatio n between the predicted ionising flux 
([Morgan et al.ll2004 ) and the internal source luminosity for 
'B' and 'C type rims is presented in FigjS] 

It would then appear that there is a link between the 
strength of the ionisation field in which the BRCs are em- 
bedded and the luminosity of stars that are produced within 
those clouds, at least for 'B' and 'C type rims. This is a 
very good indication that, at some level, the illuminating 
source of the HII region plays an important role in the for- 
mation of stars within BRCs. 



4.6. The Mass Function of tlie SFO Catalogue 

The 'Initial Mass Function' (IMF) is a key diagnostic in 
distinguishing between various regimes of star formation 
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Fig. 4. SCUBA 850/ini c ontours overlaid on DSS images of SFO 13. Triangles indicate infrared sources from the 2MASS 
Point Source Catalogue (jCutri et al.l 120031 ) that have been identified as YSOs. Contours start at 6ct and increase in 
increments of 20% of the peak flux. 
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Fig. 5. SCUBA 850/im c ontours overlaid on DSS images of SFO 12. Triangles indicate infrared sources from the 2MASS 
Point Source Catalogue (jCutri et al.l 120031 ) that have been identified as YSOs. Contours start at 4ct and increase in 
increments of 20% of the peak flux. 

Table 6. Mean core properties grouped as rim type. Values in parentheses are the standard deviation over the sample. 



Rim Type 



Tho 

(K) 



A 
B 

C 



165 (45) 
157 (29) 
127 (1) 



Tcold 

(K) 



24 (4) 
24 (4) 
22 (3) 



M 
(Ms; 



Logio (n) 
„„.-3 



Ak 



5.4 (12.5) 

11.7 (11.5) 

2.5 (1.9) 



4.8 (0.5) 
5.0 (0.4) 
5.2 (0.5) 



1.1 (1.1) 
1.5 (0.9) 

1.4 (1.2) 



processes. The Salpeter mass function (|SalpeteiJll955l ) de- (|1955D . iBlitj (jl993[ ) found an observed mass spectrum for 
scribing main sequence stars is: 'clumps' within Giant Molecular Clouds (CMC) that obeys 



dN{M) 
dM 



(X M 



-2.35 



(4) 



This relation has been refined in different circumstances 
in an attempt to describe stars that are observe d to have 
number functions differing from that presented in lSalpeted 



dN(M) 
dM 



ex Af 



-1.54 



(5) 



The masses that have been found from Eq. [2] for the de- 
tected cores in this survey are presented in Table IH the 
corresponding mass spectrum arising from this sample is 
presented in Figure [H The Salpeter and Blitz slopes are 
presented for comparison and have values of -2.35 and - 
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Fig. 6. SCUBA 850/ini c ontours overlaid on DSS images of SFO 36. Triangles indicate infrared sources from the 2MASS 
Point Source Catalogue (jCutri et al.l 120031 ) that have been identified as YSOs. Contours start at 9ct and increase in 
increments of 20% of the peak flux. 



Table 7. Comparison of cloud properties between those with and those without associated radio emission, values in 
parentheses represent the standard deviation within the group. 



IBL Detection" IBL Non-Detection 



Source Property 

THot(K) 

Tcoid(K) 

Core Mass(M0) 

Log Density 

Ak 

Nci'' 

Notts 

Lboi(L0) 

J-Jsubinm/ J-^bol 



Mean Value 


Mean Value 


154 (38) 


157 (34) 


24(4) 


23 (3) 


11 (14) 


2(2) 


5.0 (0.5) 


4.8 (0.4) 


1.5 (1.2) 


0.7 (0.5) 


3(2) 


3(3) 


1(1) 





651 (1547) 


50 (53) 


12e-03 (7.7e-03) 


14e-03 (6.4e-03) 



^Detection at 21cin CI — Class I protostcllar source 



1.54 respectively, while a least-squares fit to 
a slope of -1.67. 



our data gives well be interpreted as a lack of low-mass sources rather than 
a preponderance of high-mass sources. 



We estimate that our sample is complete to 3a at ^0.5 
Mq. Therefore, when fitting a line to the mass spectrum 
presented here we have discarded sources with derived 
masses below this value. The usefulness of a comparison 
with the Salpeter slope may not be completely valid as 
the Salpeter slope describes stars upon the main sequence, 
whereas this sample is comprised of YSOs (or possibly mul- 
tiple YSOs) which are fundamentally different astronomi- 
cal objects to main sequence stars. However, an important 
question is addressed by a comparison of the pre sented sam- 
ple's m ass spectrum with the Salpeter slope. lAlves et al.l 
(|2007l ) found that the mass function of dense molecular 
cores is similar in shape to the stellar IMF but may be 
scaled to higher masses. Their observations of 159 dense 
cores revealed a 'Salpeter-like power law'. It can be seen 
that the slope of our mass spectrum is in fact significantly 
shallower than that of the Salpeter function and somewhat 
shallower than the Blitz 'clump' function. From this, we 
can determine that these sources are more efficiently pro- 
ducing higher mass stars than may be expected or, as is 
more likely, these sources contain embedded clusters of 
stellar sources. This increased efficiency is apparent at a 
given mass, therefore, the observed difference in slope may 



The sensitivity of IRAS limits the sample to sources 
above ~ 3 L© at a mean distance of 1.29 kpc. By plot- 
ting the luminosity/distance relationship of these sourceqj 
(Fig. [ro]) it can be seen that the sample lies well above this 
sensitivity limit 



4.6.1. Separating the Sample 



If, as suggested in Section 14.51 the protostars found within 
the different cloud morphologies have some fundamental 
difference in their formation mechanism then it may become 
more evident by separating the sample on a morphological 
basis and plotting mass spectra separately. 

Here we compa r e our source sample with data taken 
from lMitchell et al.l (|200lD . in which 850 /zm SCUBA maps 
are presented of 67 discrete continuum sources in the 
Orion B molecular cloud, some of which are shown to be 



® In creating this plot only 60 and 100 /xm fluxes were used 
and compared to the flux limits of IRAS at these wavelengths. 
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Fig. 7. Histogram plot of the mass distributions of those sources both detected (filled) and not detected (black line) in 
the 21cm NVSS 
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Fig. 8. Plot of the internal luminosity of protostars vs. the predicted ionising flux arriving at the bright rim of sources 
with 'B' and 'C type morphologies 



star-forming. iMitchell et al.l (|200lh define a mass for each to the recipr ocal of the ma s s conv ersion factor d, (see 
'clump' through assuming a dust temperature, Td, via: Section 13. 2p . IMitchell et al.l (|200lh use fluxes for each 

clump taken from a 21" aperture and assume a value for 
fnK\ f Ks5o \^^ ,^§50 of 0.01 cm g^^ to find masses for a small sample of 



M, 



„p=1.50xS'85o[ea;p ( -7j^\ 



l]x 



0.01 crn^ g 



-^^od&ie clumps in their survey. In order to find masses for the 
complete sample of their observations we have recalculated 



where Sgso is the 850 ^m flux enclosed within the their dust masses assummg a dust temperature of 20 K As 
clump boundary measured in Janskys and ^550 is the mass ^^^^ 850 ^m fluxes withm a 21 aperture are not avaflable 
absorption coefficient at 850 ^lm, functionally equivalent ^^^ ^^ mamtam consistency with the method presented 
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Fig. 9. The mass spectrum of the sources detected by SCUBA in this survey. The Y-axis represents the cumulative 
number of sources (Ncum) per massbin of size 0.05 Mq. The long-dash line represents the giant molecu lar cloud clum p 
mass function taken from iBlita (jl993[ ) and the short-dash line represents the well-known Salpeter slope (jSalpeteilll955[ ). 
The solid line represents a best fit to the data from the presented sample. 



here, we have used the total integrated flux along with a 
different value of K85o=0.02 cm^g~^. The masses found 
here a gree well with those p resented for the smaller sample 
within (Mitchell et all (|200lD . 

Plo tting the mass spect rum of the clumps observed 
within [Mitchell et al.l ()2001l ) may provide a more illumi- 
nating comparison to our sample than the Salpeter func- 
tion, or the Blitz observed clump mass spectrum. As well as 
the protostellar clumps being similar objects to this sam- 
ple (the important difference being the illumination of the 
clouds by nearby massive stars) the methods used to derive 
their masses are identical, whereas the Salpeter and Blitz 
spectra are based upon various observations and assump- 
tions. It should be noted that Orion B is a single distance 
sample, while this sample represents BRCs at a range of 
distances which may not be well known. This could possi- 
bly introduce errors into the sample. As with our own mass 
spectrum lower mass sources have not been used in a best- 
fit slope to the mass spectrum in order to reduce the effects 
that an incomplete sample may have on the determination. 

The mass spectra presented in FigHT] show that clouds 
with type 'A' morphologies appear to follow a ma s s func - 
tion close to that of the survey of [Mit c hell et al.l (J200lh . 
while clouds with type 'B' and 'C morphologies have a 
significantly shallower mass function. It is noticeable that 
the best fit line to our 'A' type rim mass function may be 
skewed by a few outliers at around 3-4 Mq . The number of 
data points in this plot is low and thus we must try to fill in 
the graph around these points in order to explore this func- 
tion further. However, the suggestion from these spectra is 



that 'A' type rims seem to follow a mass law similar to 
that found in other, non-triggered regions. Whereas type 
'B' and type 'C rims have a significantly shallower mass 
spectrum, indicating a tendency toward greater numbers 
of higher mass sources. 

5. Discussion 

Our survey of the submillimetre emission found within 
BRCs has revealed the presence of embedded cores. The po- 
sitions of peak submillimetre flux are consistent with IRAS 
PSC sources, previo usly discussed by SF091, S094 and 



[Morgan et ah (I2004D and s i tuate d near the edge of known 
HII regions. [Morgan et al.l (|2004f ) established the presence 
of IBLs and photon dominated regions (PDRs) at the in- 
terface between the HII regions and t he molecular clouds 
composing the BRCs. Previous studies ([Morgan et al.ll2004 
[Thompson et al.[ [2004bt [Urguhart et all [20061 [200% also 
show that many of these clouds are likely to be in a post- 
shocked state and that any observed star formation is there- 
fore a possible result of the cloud's interaction with the 
external ionising source. 

The embedded cores are likely to contain class pro- 
tostars with few e xceptions (29 of 34 sources obey the 
[Andre et al.[ (|l993t ) classification for class protostars of 
Lsubmm/Lboi ^ 5 X 10~^). Dust tcmpcraturcs found towards 
the cores exc eed those found towards starless cores (~ 10 K 
(|Evansl[l999[ )) with an average of 24 K. This combination of 
observed submillimetre flux excess and high dust tempera- 
tures clearly establishes a scenario of ongoing early star for- 
mation within the sample of BRCs. Bolometric luminosities 
of the sources have a large range, from 1 Lq for SFO 4 to 
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Fig. 10. The luminosities of the SCUBA sample plotted as a function of distance (diamonds), the solid line represents 
the sensitivity limit of the IRAS. 
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Fig. 11. The mass spectrum of the 'A' type (left) and 'B' and 'C type (right) sources detected by SCUBA in this survey 
(pluses) and in the survey of iMitchell et al.l (|200ll ) (diamonds). The Y-axis represents the log of the cumulative number 
of sources (Ncum) per mass bin of size 0.05 M© as a fraction of N, where N is the total number of sou rces. The solid 
line is a best fit to our data above the mean sensitivity limit, the dash-dot line is a best fit to the data of IMitchell et alJ 
(|200lh . The outlier of SFO 14 with mass 51.79 M© has been neglected from the fit for the 'A' type rims as it appears 
anomalous, see Section. [5] 



6945 L0 for SFO 30. However, the majority of sources have 
Ltoi > 10 Lq suggesting that the sample consists largely 
of intermediate to high-mass star- forming regions. Though 
some, if not all, of the higher luminosity sources are likely 
host to protostellar clusters rather than individual proto- 
stars. 

Bright-rimmed clouds are the potential results of RDI 
and it has been suggested by theoretical investigations that 
they may follow a morphological evolutionary sequence 
(|Lefloch fc Lazarefj Il994l iKessel-Devnet fc Burker^ l2003t 



iMiao et aLll2006l) . In this scenario different morphological 
classifications may host different classes, or epochs, of star 
formation. A preliminary analysis revealed that this was not 
the case. The sample was separated into each morphological 
grouping and average masses, luminosities, dust tempera- 
tures, extinctions and densities were found to be similar for 
all sources with no significant deviations; no link between 
rim morphology and embedded core properties could be 
made. This strongly suggests that the development of rim 
type and the ongoing star formation within the BRCs are 
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not linked, or that they occur upon differential timescales 
such that, once initiated, star formation proceeds indepen- 
dently of any morphological changes occuring to the BRC 
as a whole. 

The detection of radio free-free emission a ssociated with 
a larg e number of the clouds in our sample ([Morgan et al.1 
I2004D establishes the presence of IBLs at the edges of many 
of these clouds. This indicates that there is some interac- 
tion occuring between the clouds and their environment, 
even if this interaction is unconnected to the internal pro- 
tostellar activity. It was found that clouds associated with 
radio-detected IBLs contain protostars, or protostellar clus- 
ters, with higher masses, densities and extinctions (column 
densities) than those not detected at radio wavelengths. 
This tendency was not found to be significantly large to 
definitively separate our sample and may simply be due to 
skewing by a few, high-mass, sources. However, when the 
predicted strength of the ionisation field within which the 
clouds are situated was investigated it was found that a sig- 
nificant correlation exists between the strength of this field 
and the luminosity of the embedded protostars or protostel- 
lar clusters. Furthermore, it was found that this correlation 
arises solely from the 'B' and 'C morphological type rims. 
This clearly indicates the interaction of the external ioni- 
sation field with star formation occuring within the clouds, 
at least in 'B' and 'C morphological type rims. The ab- 
sence of any similar correlation in 'A' type rims provides 
reasonable evidence that these clouds are not subject to 
RDI triggering processes. 

The evidence of interaction between ambient ionisation 
field strength and star formation does not rule out the possi- 
bility of formation via the 'co l lect a nd collapse' method pro- 
posed bv lElmegreen fc Ladal (|l977t) . In this scenario the ex- 
panding nebula around the star ionising the HII region gen- 
erates gravitationally unstable, shocked l ayers of gas which 
are f ragmented and have high masses (jWhitworth et all 
I1994D . This mechanism is general ly long-lived and acts upon 
large scales ("Zavagn o et al.ll2006l) and is capable of produc- 
ing the fiat edged structures seen in 'A' type rims such as 
SFO 41 and SFO 88 CFigs lXsSl fc 1X391) . 

It appears that star formation occuring within different 
morphological types of BRC follows different formation pro- 
cesses. The star formation found within 'B' and 'C morpho- 
logical rim types resulting in some part from the infiuence of 
nearby massive stars, while the formation found within 'A' 
type rims evolves presumably independently of the environ- 
ment within which the BRCs are located. This scenario con- 
tradicts the evolutionary scenario of rim type pr edicted by 
iLefloch fc Lazareff (1994) . Kessel-Devnet fc Burkert (2003 ) 
and lMiao et al.l (|2006t ). at least as connected to protostellar 
development. This inherently implies that any star forma- 
tion occuring within 'A' type rims must either have existed 
before the ionisation of the rim or that any rim ionisation 
occuring is not sufficiently dramatic to infiuence any sub- 
sequent star formation. 

In the light of the differences between 'A' type rims 
and 'B' and 'C type rims, we sought to further solidify 
any evidence of this exterior infiuence upon embedded star 
formation. An analysis of the mass functions of the two 
sub-classes of BRC revealed that the star formation within 
'A' type rims follows a similar function to protostars found 
in an environment not subjected to exterior ionisation in- 
fluence. However, 'B' and 'C type rims appear to follow 
a mass function more pre-disposed toward a higher effi- 



ciency at higher masses (see Fig fTTj) . This indicates that, in 
'B' and 'C type rims, relatively more massive objects are 
formed than low-mass stars in comparison to the Salpeter 
IMF and thus may make an important contribution to the 
overall intermediate to high-mass stellar population. 



6. Conclusions 

Our main findings are as follows: 

1. The SCUBA observations presented here generally show 
centrally condensed cores interior to the rims of BRCs. 
The morphology of the cores themselves suggests the 
presence of bound condensations, indicative of star for- 
mation. The morphology of the larger area, incorporat- 
ing the entire bright rim, is, in general, supportive of 
the scenario seen in RDI models; a dense core at the 
head of an elongated column. Support for the interac- 
tion of multiple layers is given by observations at mul- 
tiple wavelengths i.e. radio observations tracing the hot 
ionised layer bordering the FDR observed in the near in- 
frared which borders the molecular material, the dense 
regions of which are seen at submillimetre wavelengths. 

2. The physical properties of the submillimetre cores were 
determined using modified blackbody analyses. The 
dust temperatures found for the sources are high com- 
pared to starless cores (average Td = 24 K) , indicating 
the presence of internal stellar sources in these clouds. 
Based on their Lsubmm/Lboi ratio these sources are con- 
sistent with being Class protostars. In addition, the 
indicated stellar sources embedded within the clouds 
may be classified via the hypothesis that all illumina- 
tion is due to a single illuminating source. This leads to 
spectral classifications of intermediate to high mass stel- 
lar sources within the cores, whether these sources are 
isolated singular sources or weaker sources, distributed 
in clusters. 

3. Observations of mid-infrared sources taken from the 
2MASS further support the proposed scenario of early 
star formation across our sample by revealing the pres- 
ence of young stars in the vicinity of the majority 
(89%) of our sample. The spectral type of the associ- 
ated sources is not determinable from the current data, 
however the observed stars are known to be in the rel- 
atively early stages of their formation. 

4. Both the mass functions and correlations of source lumi- 
nosity with ionisation field strength suggest some fun- 
damental difference between the subsets of 'A' type rims 
and 'B' and 'C type rims. It is suggested that 'B' and 
'C type rims represent true 'triggered' star formation as 
enabled by the RDI process. There is no reason to rule 
out the possibility that 'A' type rims are also the result 
of triggering, though perhaps via a different mechanism 
i.e. th e 'collect and collapse' model (JElmegreen fc Ladal 
I1977D . The correlation between ionisation field strength 
and source luminosity is highly indicative of 'B' and 'C 
type rims forming with at least some degree of interac- 
tion with their UV surroundings. By increasing the size 
of our sample and incorporating other BRCs that have 
yet to begin their star formation this hypothesis may be 
investigated further. 
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Appendix A: Description Of Plots and Images 

Data for all of the SFO clouds are presented in the following 
pages. SCUBA 850 /xm and 450 /im contours are overlaid on 
DSS image s. Infrared sources from the 2MASS Point Source 
Catalogue (jCutri et al.ll2003[ ) that have been identified as 
YSOs are shown as triangles. 

J-H versus H-Kg colours of the 2MASS sources asso- 
ciated with each cloud are plotted along with a figure of 
the SED plot of the object composed from a best fit to 
various observed fluxes where possible. J-H versus H-Kg 
colour diagrams have solid lines representing the unred- 
dened loci of main sequence and giant stars from lKoornneeJ 
(119831). Dott ed lines represent the classical T-Tauri locus of 
iMever et al.l (1997 ) . Reddening tracks arc shown by dashed 
lines and the dash-dot line represents the reddening track 
of the earliest main-sequence stars. 
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Fig. A.l. Plots and images associated with the object SFO 1. The top images show SCUBA 450 /im (le ft) and 850 
/xm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /xm contours start at 4(T and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at Zu and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 2. Plots and images associated with the object SFO 2. The top images show SCUBA 450 /im (le ft) and 850 
/xm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003[) that have been identified as YSOs are shown as triangles. 850 /xm contours start at 9a and increase in increments 
of 20% of the peak flux value, 450 //m contours start at 5a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best flt to various observed fluxes. 
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Fig. A. 3. Plots and images associated with the object SFO 3. The top imag e shows SCUBA 8 50 /im contom'S overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 
850 /um contours start at 3a and increase in increments of 20% of the peak flux value. The bottom left plot shows the 
J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED 
plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 4. Plots and images associated with the object SFO 4. The top images show SCUBA 450 /im (le ft) and 850 
Jim, ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003() that have been identified as YSOs are shown as triangles. 850 /xm contours start at 4(T and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 5a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloudwhilc the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 5. Plots and images associated with the object SFO 5. The top images show SCUBA 450 /im (le ft) and 850 
Jim, ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003() that have been identified as YSOs are shown as triangles. 850 /xm contours start at 6a and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 3a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 6. Plots and images associated with the object SFO 7. The top images show SCUBA 450 /im (le ft) and 850 
/xm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003() that have been identified as YSOs are shown as triangles. 850 /xm contours start at 2>a and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at ia and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 7. Plots and images associated with the object SFO 9. The top images show SCUBA 450 /im (le ft) and 850 
/xm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003[) that have been identified as YSOs are shown as triangles. 850 /xm contours start at 6a and increase in increments 
of 20% of the peak flux value, 450 //m contours start at 3a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best flt to various observed fluxes. 
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Fig. A. 8. Plots and images associated with the object SFO 10. The to p image shoys 85 /im contom'S overlaid on a 
DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003[ ) are shown as triangles. 850 
/xm contours start at Scr and increase in increments of 20% of the peak flux value. The bottom left plot shows the J-H 
versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED plot 
of the object composed from a best fit to various observed fluxes. 
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Fig. A. 9. Plots and images associated with the object SFO 12. The top images show SCUBA 450 /im (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /zm contours start at 4(T and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at ia and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best flt to various observed fluxes. 
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Fig. A. 10. Plots and images associated with the object SFO 13. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 6(7 and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at ia and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 11. Plots and images associated with the object SFO 14. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
l2003f ) that have been identified as YSOs are shown as triangles. 850 /um contours start at 15it and increase in increments 
of 20% of the peak flux value, 450 ^m contours start at 3(t and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 12. Plots and images associated with the object SFO 15. The top ima ge shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 
850 /um contours start at 3a and increase in increments of 20% of the peak flux value. The bottom left plot shows the 
J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED 
plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 13. Plots and images associated with the object SFO 16. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 9(7 and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at Aa and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 14. Plots and images associated with the object SFO 18. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 9(7 and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 4(7 and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 



L.K.Morgan et al.: A SCUBA survey of bright-rimmed clouds, Online Material p 17 




QB" EJ " a^'' DO 22'" 55' 

RIGKT ASCENSIdN 
EEbtTCR: ^A. 06 2fi 58.00 DEC 423 TO T1 .7 JZOOO 



06" 23'" 05^ 00'' 22'" 55^ 

RIGHT ASCENSION 
CENTER: R A. 06 22 58,00 DEC +23 1 





Fig. A. 15. Plots and images associated with the object SFO 23. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at Su and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at ia and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 16. Plots and images associated with the object SFO 24. The left imag e shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003[ ) are shown as triangles. 
850 /im contours start at 3cr and increase in increments of 33 of the peak flux value The right plot shows the J-H versus 
H-Ks colours of the 2MASS sources associated with the cloud. 
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Fig. A. 17. Plots and images associated with the object SFO 25. The top images sho-w SCUBA 450 ^m (l eft) and 850 
Jim ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003[ ) that have been identified as YSOs are sho'wn as triangles. 850 /J,m contours start at 9(T and increase in increments 
of 20% of the peak flux value, 450 /im contours start at 2cr and increase in increments of 20% of the peak flux value. The 
central plot sho-ws the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom images 
show the SED plots of the separate core objects composed from a best fit to various observed fluxes. 
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Fig. A. 18. Plots and images associated with the object SFO 26. The top ima ge shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 
850 jjLxa. contours start at Su and increase in increments of 20% of the peak flux value. The bottom left plot shows the 
J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED 
plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 19. Plots and images associated with the object SFO 27. The top ima ge shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 
850 /um contours start at 8a and increase in increments of 20% of the peak flux value. The bottom left plot shows the 
J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED 
plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 20. Plots and images associated with the object SFO 29. The top ima ge shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 
850 jjLxa. contours start at Su and increase in increments of 20% of the peak flux value. The bottom left plot shows the 
J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED 
plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 21. Plots and images associated with the object SFO 30. The top images show SCUBA 450 //m (l eft) and 850 
yxm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003[) that have been identified as YSOs are shown as triangles. 850 ^m contours start at 12a and increase in increments 
of 20% of the peak fiux value, 450 ^m contours start at 3cr and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 22. Plots and images associated with the object SFO 31. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at Su and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at ia and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 23. Plots and images associated with the object SFO 32. The left i mage shows SCUB A 850 /zm contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003[ ) are shown as triangles. 850 
/im contours start at 3cr and increase in increments of 20% of the peak flux value. The right plot shows the J-H versus 
H-Ks colours of the 2MASS sources associated with the cloud. 
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Fig. A. 24. Plots and images associated with the object SFO 33. The top ima ge shows SCUBA 8 50 fim contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003[ ) are shown as triangles. 
850 /im contours start at ia and increase in increments of 33 of the peak flux value. The bottom left plot shows the J-H 
versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED plot 
of the object composed from a best fit to various observed fluxes. 
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Fig. A. 25. Plots and images associated with the object SFO 34. The top images show SCUBA 450 /xm (l eft) and 850 
fj,ra ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 6a and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 3a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best flt to various observed fluxes. 
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Fig. A. 26. Plots and images associated with the object SFO 35. The left imag e shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 
850 /xm contours start at 6(7 and increase in increments of 33 of the peak flux value. The right plot shows the J-H versus 
H-Kg colours of the 2MASS sources associated with the cloud. 
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Fig. A. 27. Plots and images associated with the object SFO 36. The top images show SCUBA 450 //m (l eft) and 850 
yxm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003[) that have been identified as YSOs are shown as triangles. 850 //m contours start at 9a and increase in increments 
of 20% of the peak flux value, 450 //m contours start at 9(7 and increase in increments of 20% of the peak flux value. The 
bottom plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud. 
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Fig. A. 28. Plots and images associated with the object SFO 37. The top images show SCUBA 450 /xm (l eft) and 850 
fj,ra ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 3a and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 2a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 29. Plots and images associated with the object SFO 38. The top images show SCUBA 450 /xm (l eft) and 850 
fj,m ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 9a and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 6a and increase in increments of 20% of the peak flux value. The 
bottom plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud. 
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Fig. A. 30. Plots and images associated with the object SFO 39. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
l2003f ) that have been identified as YSOs are shown as triangles. 850 /im contours start at Bit and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at Zu and increase in increments of 20% of the peak flux value. The 
central plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom images 
show the SED plots of the separate core objects composed from a best fit to various observed fluxes. 
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Fig. A. 31. Plots and images associated with the object SFO 40. The top images show SCUBA 450 /xm (l eft) and 850 
/xm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
l2003f) that have been identified as YSOs are shown as triangles. 850 /xm contours start at 4(T and increase in increments 
of 33 of the peak flux value, 450 /xm contours start at Sa and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 32. Plots and images associated with the object SFO 41. The top ima ge shows SCUBA 8 50 fim contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003[ ) are shown as triangles. 
850 /im contours start at 7a and increase in increments of 33 of the peak flux value. The bottom left plot shows the J-H 
versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED plot 
of the object composed from a best fit to various observed fluxes. 
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Fig. A. 33. Plots and images associated with the object SFO 42. The top ima ge shows SCUBA 8 50 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003[ ) are shown as triangles. 
850 /im contours start at 5(7 and increase in increments of 33 of the peak flux value. The bottom left plot shows the J-H 
versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom right image shows the SED plot 
of the object composed from a best fit to various observed fluxes. 
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Fig. A. 34. Plots and images associated with the object SFO 43. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at Su and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at ia and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 35. Plots and images associated with the object SFO 44. The top images show SCUBA 450 //m (l eft) and 850 
Jim, ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et alj 
|2003[) that have been identified as YSOs are shown as triangles. 850 //m contours start at 6a and increase in increments 
of 20% of the peak flux value, 450 //m contours start at 3a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best flt to various observed fluxes. 
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Fig. A. 36. Plots and images associated with the object SFO 45. The left i mage shows SCUB A 850 /im contours overlaid 
on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.ll2003f ) are shown as triangles. 850 
/xm contours start at 3cr and increase in increments of 20% of the peak flux value. The right plot shows the J-H versus 
H-Kg colours of the 2MASS sources associated with the cloud. 
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Fig. A. 37. Plots and images associated with the object SFO 87a. The top images show SCUBA 450 ^m (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /xm contours start at la and increase in increments 
of 20% of the peak flux value, 450 ^m contours start at 3(7 and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 38. Plots and images associated with the object SFO 87b. The top images show SCUBA 450 /im (l eft) and 850 
jUm ( right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
120031) that have been identified as YSOs are shown as triangles. 850 /xm contours start at 4(T and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 3a and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best flt to various observed fluxes. 
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Fig. A. 39. Plots and images associated with the object SFO 88. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 9(7 and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at Aa and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 
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Fig. A. 40. Plots and images associated with the object SFO 89. The top images show SCUBA 450 /xm (l eft) and 850 
um_(right) contours overlaid on a DSS image, infrared sources from the 2MASS Point Source Catalogue (jCutri et al.l 
|2003|) that have been identified as YSOs are shown as triangles. 850 /im contours start at 9(7 and increase in increments 
of 20% of the peak flux value, 450 /xm contours start at 4(7 and increase in increments of 20% of the peak flux value. The 
bottom left plot shows the J-H versus H-Kg colours of the 2MASS sources associated with the cloud while the bottom 
right image shows the SED plot of the object composed from a best fit to various observed fluxes. 



